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ABSTRACT 

Context. The term ‘solar tornadoes’ has been used to describe apparently rotating magnetic structures above the solar limb, as seen in 
high resolution images and movies from the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO). 
These often form part of the larger magnetic structure of a prominence, however the links between them remain unclear. Here we 
present plasma diagnostics on a tomado-like structure and its surroundings, seen above the limb by the Extreme-ultraviolet Imaging 
Spectrometer (EIS) aboard the Hinode satellite. 

Aims. We aim to extend our view of the velocity patterns seen in tornado-like structures with EIS to a wider range of tempera¬ 
tures and to use density diagnostics, non-thermal line widths, and differential emission measures to provide insight into the physical 
characteristics of the plasma. 

Methods. Using Gaussian fitting to fit and de-blend the spectral lines seen by EIS, we calculated line-of-sight velocities and non- 
thermal line widths. Along with information from the CHIANTI database, we used line intensity ratios to calculate electron densities 
at each pixel. Using a regularised inversion code we also calculated the differential emission measure (DEM) at different locations in 
the prominence. 

Results. The split Doppler-shift pattern is found to be visible down to a temperature of around log T = 6.0. At temperatures lower 
than this, the pattern is unclear in this data set. We obtain an electron density of logne = 8.5 when looking towards the centre of 
the tornado structure at a plasma temperature of log T = 6.2, as compared to the surroundings of the tornado structure where we 
find logne to be nearer 9. Non-thermal line widths show broader profiles at the tornado location when compared to the surrounding 
corona. We discuss the differential emission measure in both the tornado and the prominence body, which suggests that there is more 
contribution in the tornado at temperatures below log T = 6.0 than in the prominence. 
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1. Introduction 


The determination of plasma properties is an essential com¬ 
ponent of our understanding of the structures observed in the 
solar atmosphere, and it provides important constraints on the 
scenarios attempting to explain their properties and appearance 
( [Labrosse et al.|2QI^ . This work aims to provide an insight into 
the physical conditions found in the so-called large-scale solar 
tornadoes that have been observed in the legs of prominences 
( [Wedemeyer et al.|2QI3| ). There is an ongoing debate about how 
to interpret the apparent rotation of these features (see [Orozco 


Suarez et al.||2QI2| O et al.||2QI2| |Su et al.||20I2| [Panesar 


2QI3[ |Su et al.||2QI4| ). Clearly, combining high-resolution im¬ 
ages, spectroscopic observations, and magnetic field measure¬ 
ments is the best approach to making any progress in this dis¬ 
cussion. 

Unfortunately, direct measurements of the magnetic field 
in tornado-like prominences are difficult to perform. It is log¬ 
ical then to combine imaging and spectroscopic data analysis 
to unde rstand the nature of these a p parent motions. This wa s 
Orozco Suarez et al.] ( |20 1 2| ), [Wedemeyer et al.| ( [2QI3| ), 


done by 


Su et al 


( |20I4| ), amongst others, and their results suggest that 


the observed prominence foot structure is rotating. However, 


their analysis leaves open a few questions related to the state 
of the plasma that was observed. In particular, the temperature 
and density of the plasma are two important pieces of informa¬ 
tion that are needed to put forward and test physical models of 
tomado-like prominence stmctures. Similarly, any signature of 
non-thermal processes in the structure may give a clue as to what 
types of mechanisms are at work. 

In this paper, we investigate the plasma properties of a so¬ 
lar tornado using data obtained by Hi node/EIS and SD O/AIA on 
14 September 2013 (also studied by |Su et al.|20I4| ). Section 
gives an overview of the observations and explains the proce¬ 
dure followed for the data analysis. The analysis of the spectral 
lines used in this work allowed us to make line-of-sight velocity 
measurements (Section |^, density diagnostics (Section]^, and 
a determination of non-thermal line widths in the observed re¬ 
gion (Section [^. We present our differential emission measure 
analysis in Section Section]^ gives our conclusions. 


2. Observational overview 


The Extrem e-ultraviolet Imaging Spectrometer (EIS [Culhane 


et al.| (2007] ) aboard the Hinode spacecraft uses a slit (I” or 
width) or slot (40" or 266" width) followed by a grating 
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spectrometer to pass EUV light to two CCDs (long waveband 
- 246-292 A - and short waveband - 170-211 A). A few studies 
have been used to observe prominences with EIS, but very few 
have focussed on tornadoes. Two related EIS observing plans 
named eis_tornadoes_scan and eis_tornadoes_sns were 
designed by P. Gomory et al. to look for tornadoes above the 
limb, and they consisted of a raster and a sit-and-stare observa¬ 
tion, respectively. The study made use of the 2" slit of EIS, tak¬ 
ing 50 slit positions and covering a spatial extent in x of 100". 
The y extent of the raster was 256". 


The study was run over seven days in September 2013 (9, 11, 
12, 13, 14, 18 and 19, beginning at around 03:30UT each day) 
with each day’s observing plan consisting of a raster followed 
by a three hour sit-and-stare, and hnahy a second raster. One of 
the seven day’s observations (14 Sept.) caught a tornado with 
the hrst raster and subsequent sit-and-stare. It is this data set that 
comprises the focus of this paper. By the second raster of the 
day, the tornado was no longer visible to EIS. 


This data set was also used by |Su et ST] ( |2014| ). In that pa¬ 
per the authors discussed the line-of-sight velocities observed in 
coronal lines formed at around 1.5-2 million K, concluding that 
the hot plasma observed in both the raster and sit-and-stare indi¬ 
cates that the tornado structure is rotating. 


2.1. Description of the event 

Eigure shows the tornado as seen in the 171 A waveband of 
AIA (central panel), as well as Solar Magnetic Activity Re¬ 
search Telescope (SMART, Hida Observatory, Kyoto University, 
Japan) H-a (left panel) and the larger prominence structure in 
AIA 304 A (right panel). In coronal lines the tornado is visi¬ 
ble to a height of around 20" as a dark feature above the limb. 
It appears darker than the surrounding corona because there is 
a certain volume occupied by cool plasma that is not emitting 
in hot coronal lines (emissivity blocking), and because the cool 


plasma abs orbs the radiation coming from behind (see [Heinzel 
|et al.|200^ and references therein). We know that the dark col 
umn on the limb contains hydrogen and helium. The presence of 
H I, He I, and He ii will lead to absorption of coronal radia tion at 


wavelengths below 912 A, 504 A, and 228 A respectively (Orrall 
|& Schmahl|1976[[Labrosse et al.|2010| ). In these dark structures, 
nearly ah background emission in coronal lines in the EIS spec¬ 
trum is therefore absorbed by the cool plasma along the line of 
sight ( [Labrosse et al.|20lT] ). Therefore, the emission observed by 
EIS in the tornado comes from a sheath of hot plasma in front 
of the tornado and from the corona in front of it. Eurthermore, 
the H-a tornado is similar in size to the dark feature seen in AIA 
171 A images, which is explained by the fact that the optical 
thickness of the plasma at 171 A is comparable to that of H-a 
( [Anzer & HeinzeH2005| ). The prominence can be seen as a hl- 
ament on disc for a number of days before it crosses the limb 
(Eigurej^. 

The EIS raster has a smaller held of view than the AIA and 
SMART images (Eigure[T]). Eigure shows the intensity image 
resulting from the Gaussian htting of the 195 A Ee xii line along 
with an AIA image of the 193 A waveband. The emission from 
this AIA window comes predominantly from Ee xii. 

The 195.12 A line is the line in this EIS study that most 
clearly shows the tornado structure above the limb. In addition to 
the tornado visible above the limb, this image reveals the pres¬ 
ence of a structure extending southwards from the top of the tor¬ 
nado similar to smoke from a chimney, which is co-spatial with 
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Table 1. EIS spectral lines used in the analysis 


Ion 

Wavelength (A) 

logT(K) 

Blend 

Eex 

184.537 

6.1 

Ee XI 

Ee VIII 

185.213 

5.7 

- 

Eexi 

188.216 

6.2 

Ee XII, Ee xi 

Ee XI 

188.299 

6.2 

Ee XII, Ee xi 

Ee XI 

192.627 

6.2 

- 

Eexi 

192.814 

6.2 

0 V, Ca XVII 

Ov 

192.904 

5.4 

0 V, Ee XI, Ca xvii 

Ee XII 

195.119 

6.2 

Ee XII 

Ee XII 

195.179 

6.2 

Ee XII 

Ee IX 

197.862 

6.0 

- 

Ee XIII 

202.044 

6.3 

- 

Si VII 

275.361 

5.8 

- 


the main prominence body seen in 304 A and with the faint emis¬ 
sion seen in H-cr. 


2.2. Data reduction 


EIS data is aquired as level-0 data hies, so it is necessary to run 
the standard eis_prep procedure hrst to bring the data up to 
level-1. This procedure removes any hot or warm pixels, dark 
currents, and cosmic ray hits and performs an absolute calibra¬ 
tion of the data. We used CHIANTI ( |Dere et ar]|1997t |Landi| 
|et al.|2012| ) to identify ah the lines visible in the calibrated spec¬ 
tral windows, and ht them assuming Gaussian prohles with the 
help of Craig Markwardt’s mpfit package. The wavelength cal¬ 
ibration was carried out using standard routines to correct for 
orbital variations. We then corrected for solar rotation, choosing 
the central part of the tornado to be hxed with a line-of-sight 
velocity of zero. 


It has been noted that the detectors aboard the EIS instrument 
are decaying as a function of b oth time and wavelength ( |Del| 
|Zanna|2QT3 [Warren et al.|2014| ) instead of simply as a function 
of time, so alternative, corrective intensity calibrations must be 
employed. Two similar methods are well described in tho se pa¬ 
pers, and it is the method outlined by [Warren et al.| ( [2QT^ that is 
used here. 


2.3. EIS spectral lines 

All lines used in this analysis are listed in TableThe emission 
of these lines is from relatively hot plasma (where logT > 5.4) 
when compared to the cool plasma of the prominence (log T ~ 
4). Here we primarily work on the emission spectra observed 
when looking along the line of sight towards the tornado. Ac¬ 
cording to [Parenti et al. [ ( [TO 1 2[ ) this emission is mostly from the 
prominence-to-corona transition region (PCTR). The word ‘tor¬ 
nado’ here refers to all the hot plasma in the temperature range 
logT = 5.4 - 6.3 along this line of sight 

Due to the finite spectral resolution of the EIS instrument, 
we observe a number of blended lines. It is possible to resolve 
some of these blends. The process is described below. 

2.3.1. Fex 184.537 A 

The Eex line at 184.537 A has a small component of Eexi in 
its blue wing. Although much smaller than the Eex line itself, 
this Ee XI line is non-negligible. Due to the fact that Ee x domi- 
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SMART H-alpha 14-Sep-2013 02:32:14 UT SDO AIA_3 171 14-Sep-2013 03:22:23 UT SDO AIA_4 304 14-Sep-2013 03:22:31 UT 





Fig. 1. Images from the Solar Magnetic Activity Research Telescope (SMART, Hida Observatory) and the Atmospheric Imaging Assembly (AIA) 
aboard the Solar Dynamics Observatory (SDO) showing K-a from SMART (left), the 171 A passband (centre), and the 304 A passband (right) 
from AIA. The central image shows a tornado, absorbing background coronal plasma emission, with the right hand image giving context as to the 
larger prominence structure. Also shown in white in the central image is the field of view for the EIS raster. 
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Fig. 2. Images from SMART K-a (left), AIA 171 A passband (centre), and AIA 304 A passband (right). 


nates the blend, however, a simple two Gaussian fit successfully 
removed the Fe xi contribution from the wing. 


2.3.2. The Fexi 188.2 A blend 

The Fexi doublet at 188.2 A is a pair of lines that are partly re¬ 
solvable, making de-blending them fairly simple, however there 
is a small component of Fexii in the blue wing of the profile 
that is not negligible. A triple Gaussian does a good job of fitting 
this group, but the fit can be improved by tying the centroids of 
the Fe xi lines together. We cannot tie the intensities together as 
there is a slight density sensitivity between this pair. 


2.3.3. The Fe xi and O v 192 A blend 


This blend at 192 A represents one of the most complex blends 
seen by EIS. A number of authors have worked on de-blending 
these lines (see, e.g. |Ko et al.|2009|[Graham et al.|2013| ), mostly 


in flaring or active region rasters where the Ca xvii line (forma¬ 
tion temperature logT = 6.8) will dominate. In this raster, how¬ 
ever, we would expect a minimal amount of emission at this tem¬ 
perature. This changes our approach to de-blending these lines, 
as the Caxvii line will contribute a negligible amount to the 
blend. Instead, the Fexi line becomes the dominant line, with 
a fairly large amount of O v emission in this region too. 

In this spectral window we also have the Fexi line at 


192.627 A. This line is resolved, outside of the main blend, and 
can be fitted with a single Gaussian. Although it is weaker than 
the Fe xi line within the blend, it is always visible in this raster, 
and can be used to constrain the centroid of the blend ed com po- 
nent. As with the 188.2 A lines discussed in Section 


2.3.2 


the 

intensities of these two Fe xi lines cannot be tied together as they 
display a density sensitivity. 

The five O v lines in the blend are difficult to deal with. These 
are the only O v lines available in this study, so we have no way 
of tying parameters to other O v lines outside of the blend. These 
lines were handled by making two assumptions about them. The 
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SDO AIA_2 193 14-Sep-2013 03:22:18 UT 
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Fig. 3. Intensity maps of the AIA 193 A waveband (left) and the EIS Fexii 195.12 A spectral line obtained by Gaussian fitting (right), showing the 
tornado and the structure extending southwards above the limb. The plus markers in the EIS image (right) indicate the points at which the tornado, 
prominence and coronal DEMs presented in Sectionj^are calculated. The AIA image is using the same waveband as the online movie, that shows 
the field of view of these images as a white rectangle. 


first is that the five lines can be modelled by two Gaussians 
centered about the two strongest of the lines - 192.904 A and 
192.797 A. The second is that these two Gaussians are tied in 
both centroid position and in intensity. This is not entirely ac¬ 
curate, as these two lines normally form a density sensitive pair, 
however, at the densities considered here this sensitivity is min¬ 
imal. As in IGraham et al.| ( |2Q13| ) we here assume a fixed den¬ 
sity of 3 X 10^^ cm“^, and therefore a fixed intensity ratio of 
A 92.797 = 0. 39 / 192 . 904 , for these lines. With the Fexi line already 
removed from the blend this is a good approximation for these 
lines, even if the assumed density is higher than might be ex¬ 
pected for this line - [Labrosse et al.| ( |201Q| table 4 and refer¬ 
ences therein) suggest that the density of O v would be around 
log He ~ 9.5 for an active region prominence, slightly lower for 
quiescent prominences - the density sensitivity between these 
two lines at these densities is negligible, and the resulting fit is 
sufficient for the analysis performed here. 

Figure]^ shows the results of this de-blending process. Seen 
on the left hand side is the Fexi, 192.627 A line, shown fitted 
with an orange Gaussian. The fit for the other Fe xi line is plotted 
in green, with the two O v lines in purple and blue. The thin dark 
blue line represents the total contribution from both O v lines. 

2.3.4. The Fexii 195.1 A blend 

The Fexii 195.119 A line is the strongest fine in the EIS spec¬ 
trum, and is one of the best for analysis of this tornado. In the 
red wing of this line is a small component of Fexii at 195.179 A 
that is non-negligible. The 195.119 A line was one of those used 


Su et al.| ( |2014| ) in their analysis of this event, and here we 


by _ 

follow a similar approach to remove the red-wing line. We take 
a double Gaussian fit, and as with the Fe xi doublet in Section 


2.3.2| we tie the centroid locations together. Here again we can¬ 
not tie the intensities together to form a constant ratio as these 
two fines form another density sensitive pair. The fitted 195 A 
spectral window is shown in Figure 

3. Line-of-sight velocity measurements 

After all lines visible in the spectral windows were fitted, we 
used line centroids to calculate the lin e-of-sight veloc ity at each 
pixel. This is the same method used by |Su et ar] ( |2014| ), who sug¬ 
gest that they found evidence of rotation of this structure in lines 
formed at coronal plasma temperatures, specifically in Fe xii and 
XIII lines formed at ~ 1.5 MK and ~2 MK respectively. It is found, 
when looking at other lines available in the raster, that this tem¬ 
perature range in which we see this split pattern of redshifts and 
blueshifts along the axis of the structure can be extended down 
to at least 10^ K. 

Figurej^shows line-of-sight velocities across the tornado for 
iron lines (Feix - Fexiii), with error bars deriving directly from 
the uncertainties in the fit parameters. The fines used in this plot 
are all formed at temperatures in the range log T = 6.0-6.3. Here 
we see a Doppler split about the axis of the tornado in each of 
these lines, with the north side (left half of Figurered shifted 
and the south (right half of Figure side blue shifted. 

The Feix line at 197.86 A is formed at ~10^K, and although 
the pattern is less clear, we see in Figure [7] the same red/blue 
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Close up tornado Fe IX 197.862 



800 805 810 815 820 825 830 835 

X (arcsecs) 

Intensity (erg cm'^ s"' sr"*) 


7 17 27 38 48 58 68 


Close up tornado Fe IX 197.862 



800 805 810 815 820 825 830 835 

X (arcsecs) 

Velocity (km s'^) 


-15 -10 -5 0 5 10 15 


Fig. 7. Intensity (left panel) and line-of-sight velocity (right panel) maps for the Feix 197.86 A line. The contour plot on the velocity map is the 
Fe XII195 A line intensity (Figurej^, that shows the tornado structure most clearly above the limb. 
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Fig. 4.192 A blend, the most complex blend under consideration in this 
EIS raster. Shown here is the de-blended line, fitted using four Gaus¬ 
sian components, which are the coloured lines. The original spectrum 
is shown here as the histogram with error bars. The orange line is the 
Fexi 192.627 A line, green is Fexi 192.814 A, and the purple and blue 
lines represent the O v lines. The thin, dark blue line shows the overall 
contribution from the O v lines, and the solid black line is the total fit 
to the spectrum. The log scale on the y-axis enhances the weaker lines 
relative to the stronger lines. 


Fig. 5.195 A spectral window as seen by EIS, fitted with multiple Gaus- 
sians. Here the histogram with error bars is the raw data as measured by 
EIS, the dark red curve is the Fexii 195.119 A line and the magenta 
curve is the Fexii 195.179 A component. There is a third Gaussian fit¬ 
ted here, in red, at around 195.4 A, that is an unidentified line. The solid 
black curve is the overall profile created from these Gaussians. The log 
scale on the y-axis emphasises the fits of the weaker lines. 


pattern that is described in the |Su et ^ paper for hotter plasma 
temperatures. 

We cannot see clear evidence of this Doppler pattern for 
other lines in the study that are formed at temperatures less than 
~1 MK. For example the Fe viii line at 185.21 A, formation tem¬ 
perature of logT = 5.7 and shown in Figureshows no signs 
of the patterns visible at higher plasma temperatures. This can¬ 


not, however, be taken as evidence that there is no split Doppler 
pattern at these temperatures, as all lines available in this tem¬ 
perature range are affected by blends and/or low signal-to-noise 
ratio (SNR). For example, the O v lines are heavily blended, (see 
Section |?.3.3| ). In order to see any systematic line shifts we need 
these to be accurately de-blended, as the typical line-of-sight ve¬ 
locities are small (on the order of a few km s“^). Low SNR for 
some lines above the limb also make it difficult to measure line- 
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Velocity vs. position across tornado 



Fig. 6. Line-of-sight velocities for Fe ix - Fe xiii ions for a cut through 
the tornado - ‘cut 3’ as seen in Figure {U] The tornado has a spatial 
extent from approximately [-5,5] arcsec on the x-axis of this plot. 


Close up tornado Fe VIIM 85.213 



X (arcsecs) 
Velocity (km s'^) 



Fig. 8. Line-of-sight velocity map for the Feviii 185.21 A line. The 
contour plot is the Fexii 195 A line intensity (Figure [^. 


CHIANTI V. 7.1.3 Fe XII 195.119 (A)/1 95.179 (A) T= 1.58e+06(K) 



Fig. 9. Electron density versus intensity ratio for the Fe xii 
195.119 A/195.179 A line pair. 


three to four pixels from the intensity maximum (EIS Wiki: 
http://solarb.mssl.ucl.ac.uk: 8080/eiswiki/Wiki.jsp?page=Spatial 
%20offset%20of%20intensity%20and%20velocity%20features). 
When there is an increasing intensity gradient from north to 
south, a blueshift is introduced, whereas decreasing intensity 
gradients north to south introduce redshifts. In either direction it 
is found that the additional Doppler velocity is around 5 km s“^ 

If this instrumental line shift is present in this data set, then 
the observations presented here could be seen as spurious. How¬ 
ever, we see in the data presented in Figure |7j especially in the 
top half of the tornado, that there is not a large intensity gradient, 
yet the Doppler pattern persists. Also, as can be seen in Figure 
in the Fe x and Fe xi lines we recover Doppler shifts exceeding + 
5 km s“^ so above the level of the estimated instrumental shift. 

4. Density diagnostics 

There are a number of density sensitive lines visible to EIS, and 
a few of those were available in this study. Using the CHIANTI 
package ( |Dere et al.|1997t [Landi et al.|2012| ), that has informa¬ 
tion on expected intensity ratios for density sensitive line pairs, 
we can calculate the electron density at each pixel position. Fig¬ 
ure shows an example density curve for the Fe xii doublet 
195.119 A and 195.179 A. 

There are four available density diagnostic pairs in this study, 
three from Fexi and Fexii, all corresponding to temperatures 
around log T = 6.2. A potential O v diagnostic was ruled out be¬ 
cause of the difficulties in fitting the blended lines, as described 
in Section l233l 


of-sight velocities accurately. This is especially evident in the 
Si VII 275.36 A line (formation temperature of log T = 5.8) and 
in the Fe viii 185.21 A line. 


One effect that we must acknowledge, that could be re¬ 
sponsible for the observations of the split Doppler pattern at 
higher temperatures, is an artificial Doppler shift that may 
be introduced by a point spread function (PSF) that is both 
tilted and elliptical, similar to that seen in SOHO/CDS data 
( [Haugan|1999| ). According to [Young et al.| ( [20T^ Appendix B) 
additional Doppler shifts may be found in EIS spectra in regions 
where there is a large intensity gradient north-south along the 
slit. This instrumental line shift is typically offset by about 
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4.1. Fexi diagnostics 

The three available Fexi diagnostics all make use of the 
192.627 A line as one of the diagnostic pair. The other three 
lines used in the diagnostics are each of the 188 A pair and the 
192.814 A line. There are no useful mutual diagnostics between 
the three of these, meaning we only have three diagnostic pairs in 
total for this ion. In each of these, however, we see a slight dip in 
density at the tornado location. This is seen in Figure which 
shows one of these Fexi diagnostics with the Fexii 195.119 A 
intensity contour plotted on top. The other two density maps 
provide a similar result, showing a lower density region when 
looking towards the tornado centre. 
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Close up tornado Fe XI 192.627/192.814 
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Fe XII Density vs. position across tornado axis 



Fig. 12. Density plotted as a function of position across the tornado 
axis for five parallel cuts. Cut 1 is the uppermost as seen in Figure {Tl] 
(farthest from the limb) and cut 5 is closest to the limb. Zero on the x- 
axis is defined as the centre of the tornado, where the Doppler shift is 
nearest to zero. 


Fig. 10. Density map for one of the Fexi diagnostics. The other two 
available provide a very similar result. These are not shown here, but 
are consistent with the diagnostic presented. 


4.2. Fexii diagnostic 

The Fexii density diagnostic available consists of the doublet 
at 195 A. This is a blend that is dominated by the 195.119 A 
line, with the weaker line centred at 195.179 A contributing a 
non-negligible amount to the red wing of the blend. Using the 
fitting techniques outlined in Section |2.3.4| these lines can be 


de-blended and used for the diagnostic. 

Figure shows the resulting density map. Just like in the 
Fe XI maps, we see a small dip in density at the tornado location. 
This is seen more clearly if we take cuts through the tornado axis, 
shown as white parallel lines in the right hand panel of Figure 
the results of which are plotted in Figure 

Figureshows how the density changes across the tornado 
axis. Here we can clearly see that the density is lower towards the 
middle of the tornado, especially looking at cut 3 (triangle mark¬ 
ers, magenta line). This cut reveals a drop in density of around a 
factor of three when looking towards the centre of the tornado. A 
number of the other cuts also show a drop in density when look¬ 
ing towards the centre of the tornado, with a minimum at around 
this position. The uppermost cut (plus markers, light green line) 
does not display the same pattern, but this location is further into 
the corona and it is to be expected that the ambient electron den¬ 
sity is lower here than nearer to the limb. 


4.3. Density analysis 

From Figure p^ we observe a mean electron density of log^^ ~ 
8.5 - 9.2 in the hot (logT = 6.2) plasma. This range is consis¬ 
tent with values found previously (see review by|Labrosse et al. 

|MT0| ). 

These observations bring up questions about the origin of the 
observed drop in density. The velocity maps presented in Figures 
and [^suggest that the majority of the hot line emission is com¬ 
ing from the tornado - the strong Doppler structure indicates that 


these emission lines are formed in the moving plasma. Hence we 
can assume that this dip in density is also due to the tornado it¬ 
self. We also know that the temperature of this lower density 
region is around log T = 6.2, as it is visible in all available di¬ 
agnostics at this temperature. This suggests that the rotating hot 
plasma is less dense in the central parts of the tornado than at its 
edges, while the edges are denser than the surrounding corona. 


5. Non-thermal line widths 


The measured spectral line width can be attributed to three 
main factors. The first is the instrumental width, broadening 
introduced by the optics of the telescope, that has a measured 
value and can easily be removed. Secondly there is the thermal 
width, which is associated with the thermal motion due to ambi¬ 
ent plasma temperatures. Any remaining line width is the non- 
thermal line width (NTLW). Assuming that each of these can be 
characterised by a Gaussian profile, we can split the measured 
line width, AX^eas^ by equation 

^^meas “ ^Anst ^ ^ ^At ’ W 


where AXinst is the instrumental line width, AXth is the ther¬ 
mal width and AX^t is the non-thermal component. The thermal 
width can be calculated, and the instrumental width of EIS has 
been measured (see [Brown et al.|2()0^ . Here we adopt the mea¬ 
sured laboratory value (0.047 A) for the short wave band CCD 
as an absolute lower limit for AXinst, so that we can ascribe an 
upper limit to the non-thermal width. 

With the instrumental and thermal line widths removed from 
the line profile, we can examine the non-thermal broadening. 
Figure shows the NTLW maps of the Fexii 195.119 A, 
which was de-blended as described in Section |2.3.4[ and Fe xiii 
202.044 A lines (other lines did not yield clear NTLW maps). Er¬ 
rors were calculated using the errors on the measured line width 
from the fitting. The errors on the thermal (calculated using CHI¬ 
ANTI values) and instrumental line widths are assumed to be 
zero. 

In both cases we see a slightly broadened profile at the tor¬ 
nado location when compared with the corona immediately next 
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Fe XII 195.119/195.179 14-Sep-2013 03:02:34 UT 
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Fig. 11. Density map for the Fe xii line pair at 195 A. The right panel is a zoomed-in image, the dimensions of which are shown as a white box in 
the left hand panel. Overplotted on both is the contour of the 195.119 A line, showing the position of the tornado. Also plotted on the right hand 
image are five parallel cuts through the tornado axis, used in Figure 


Fe XII 195.119 NTLW vs. position across tornado axis 
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Fig. 14. Non-thermal line width plotted as a function of position across 
the tornado axis for five parallel cuts through the tornado in the NTLW 
map, here shown for the 195A Fexii line. These are the same cuts used 
in Figure and are shown in Figure [TT] 


to it, as is emphasised in Figure which shows cuts through 
the tornado axis in NTLW. Here we clearly see broader profiles 
when looking at the tornado. This would suggest that there is an 
additional broadening mechanism here. Although it seems that 
there is broadening local to the tornado, we cannot take this as 
conclusive evidence that the tornado is the cause of the broad¬ 
ening. As can be seen in Figure p3] there are similar patterns of 
broadening at other limb locations, so we cannot rule out that 
this is a chance occurrence. 
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If we assume, however, that the broadening is indeed caused 
by the tornado structure, we must consider what mechanisms 
could be causing this non-thermal broadening. One possibility 
is related to the possible existence of two types of magneto¬ 
hydrodynamic systems in the prominence structure: the dens¬ 
est parts of the prominence plasma would be supported by hori¬ 
zontal and relatively weak magnetic fields, while more dynamic 
parts of the prominence could be related to strong horizontal 
fields combined with a turbulent field ( [Schmieder et al.|[2014| ). 
This could explain non-thermal line broadening in the tornado. 
Alternatively, it has been suggested (T. Zaqarashvili, private 
communication) that the Kelvin-Helmholtz instability in torna¬ 
does could cause line broadening. This kind of instability occurs 
when there is a shear fiow between two plasmas, much like we 
would expect at the boundary between a rotating tornado and 
the surrounding corona. This observation of non-thermal broad¬ 
ening then appears to add weight to the suggestion that the tor¬ 
nado structure is rotating. Observations of the Kelvin-Helmholtz 
instability have been suggested in relation to observed Rayleigh- 
Taylor instabilities in quiescent prominences ( |Berger|2014| ). 

6. Emission measure distribution 

In this section we investigate the emission measure distribution 
and differential emission measure at four locations in the raster: 
in the tornado, in the main prominence body, and at two locations 
in the corona at the same altitude as the tornado and prominence 
respectively. 

Differential emission measures (DEMs) are a way of iden¬ 
tifying the temperature distribution along the line of sight. Us¬ 
ing a regularised inversion code from [Hannah & Kontar | ( |20 1 2| ) , 
adapted for use of spectral lines available in this EIS data set, 
we can calculate the regularised DEM at different points across 
the raster. This code utilises the intensities from the Gaussian 
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Fig. 13. Non-thermal line width maps of the Fe xii 195.119 A (left) and Fe xiii 202.044 A (right) lines. These images show a slightly larger field of 
view than, for example, Figures [^and|^ but the resolution is the same. 


Table 2. List of lines available in the FIS study for constraining the 
DFM. 


Ion 

^o(A) 

logT(K) 

Ov 

192.904 

5.4 

Ee VIII 

185.213 

5.7 

Si VII 

275.361 

5.8 

Ee IX 

197.862 

6.0 

Eex 

184.537 

6.1 

Eexi 

188.216 

6.2 

Ee XI 

188.299 

6.2 

Ee XI 

192.814 

6.2 

Ee XII 

195.119 

6.2 


fitting procedure described in Section [2 along with contribution 


functions c alculated using CHIANTI ( |Dere et al. 
|et al.|2Q12| ) to perform the inversion (see Hannah & 


T997llLa5di| 


■Contar|20T^ 


for full description of inversion technique). We assumed photo- 
spheric abundances for calculating the DEM in pixels located in 
the tornado and prominence body, and coronal abundances for 
pixels located in the corona. 


DEMs in prominences have been calculated numerous times 
using a number of different instruments (see, e.g. Schmahl & Or- 
rall||1986t |Wiik et al.||1993| [Parenti & Vial||2007[ |Parenti et al. 

2012| ), but none have used EIS data for this analysis. To our 
knowledge, no DEMs of solar tornadoes have been published 
at the time of writing. 

Table |2] contains information on the EIS lines that have been 
used when calculating these DEMs. Some lines have been omit¬ 
ted from the DEM for a number of reasons. Eor example, the 


coolest line in the study, Hen at 256 A, must be ignored due to 
the fact that it is optically thick, as well as being heavily blended. 
Ee XIII 202.044 A, the hottest line available in the raster, was also 
not used as it is found to be density sensitive at densities below 
around log^^ = 10 (see figure 4.2 in |Graham|2014[ ). In total, we 
are left with nine lines that can be used for the DEM. 


Eor this analysis errors of 22% have been assumed for the 
line intensities. This is considerably higher than the errors from 
line fitting, but due to combined uncertainties from the response 
function calibration of the EIS instrument ( |Lang et al.|2006| ) we 
assume here the larger percentage error. This calibration uncer¬ 
tainty was measured pre-flight, and so must be considered as a 
lower limit. 


Eigure (left panel) shows a DEM produced from a pixel 
in the tornado, and Eigure [T^ (left panel) shows that from a pixel 
in the main prominence body (indicated by plus signs in Eig.[^. 
We recover a different solution for each of these regions, with 
the tornado (Eigure[^ showing more plasma at log T = 5.4 than 
in the prominence: the DEM is larger in the tornado than in the 
prominence by a factor of 4 - 5 at this temperature. This suggests 
that there is more cool material along the line of sight towards 
the tornado, which is supported by the fact that when observed 
in coronal lines, the tornado appears darker than the prominence, 
signalling greater absorption of the background emission. 

Eigures [17] and show DEM and EMD plots for pixels in 
the corona, at similar altitudes to the DEMs presented in Eigures 
[T5| and \T6\ respectively. We here present these coronal DEMs to 
show that it is not just the foreground corona that we are seeing 
in the prominence and tornado DEMs. Eigure [T7| (for a pixel in 
the corona, close to the tornado) contains one main peak, with 
DEM(T) = 4 X 10^^ cm“^ K“\ at around logT = 6.2, and a 
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logio Temperature [K] 


Fig. 15. Result of the DEM code for a pixel in the tornado (see Figure [^. Photospheric abundances were assumed for this DEM. Plotted here 
is the resulting differential emission measure (left panel) along with the emission measure distribution (right panel). Also plotted on the EMD 
panel (right) are the EM loci curves, or constraint curves, for each of the lines used in the DEM calculation (see Tablej^for line identification and 
formation temperatures). The green ‘error bars’ in the DEM (left panel) and the grey ranges in the EMD (ri ght panel) represent con fidence regions 
for the DEM and EMD fits respectively. Plotted as pink asterisk markers on the DEM are the results from |Parenti & ViaT| ( |2007| ) for a quiescent 
prominence. 




log^o Temperature [K] 


Fig. 16. Same as Fig. 15 for a pixel in the main prominence body. 


smaller (DEM(T) ~ 1.5 x 10^^ cm“^ K“^), less well constrained 
contribution at around log T = 5.5. In contrast, the tornado DEM 
(Eigure[^ has a lower temperature peak at logT = 5.4, where 
DEM(T) = 8x10^^, suggesting that there is much more contri¬ 
bution at lower temperatures coming from the tornado itself. In 
Eigure (for a pixel in the corona at the same altitude as the 
prominence body) we see little contribution below logT = 6.0, 
consistent with the findings of previous authors ( e.g. |Landi &| 
|Young||MT0| ). The coronal DEMs (Eigures \rf \ and 1 1 8|) confirm 
that the peaks around log T = 6.2 observed in the DEMs of the 
tornado and the prominence (Eigures p~5] and could be due to 
foreground emission. 

Taking an estimate of the gradient of the slope of these DEM 
curves between logT = 5.65 and logT = 6.10 for the tornado 

Article number, page 10 of 


and logT = 5.88 and logT = 6.10 for the prominence body, we 
recover a value of 2.6^^^ for the tornado (Eigure [T^ and 6.4^^-^ 
for the prominence body (Eigure p^. 

The coloured lines in the right panels of Eigures [T^ [T7] 
and represent the constraint curves for the EMDs. Even 
though we have a limited number of lines, there is still quite 
a good temperature range available to constrain the curves. 
More lines formed at cooler temperatures would have helped, 
as prominences generally consist of lower temperature plasma. 


Plotted as asterisks in the DEM panels of Eigures p3] and \TE\ 
are the results of |Parenti & Vial| ( |2007| ) for a quiescent promi¬ 
nence. In that paper they used SOHO/SUMER data on a large 
prominence, finding a high temperature peak in the DEM at 
around logT = 6.0, by using a selection of lines from the 
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Fig. 17. DEM (left) and EMD (right) for a point in the corona at a similar altitude as the tornado DEM (Eigurep"^. 



Iogio Temperature [K] 


Fig. 18. DEM (left) and EMD (right) for an arbitrary point in the corona 

M 



a similar height above the limb to the prominence body DEM of Eigure 


SUMER spectrum. The analysis done here is on a much smaller 
prominence, only 20" in altitude, using a smaller number of 
EIS lines, most of which are formed at higher temperatures than 
those used by [Parenti & Vial| Those authors also used an older 
version of the CHIANTI database for their analysis (version 4.2, 
as opposed to version 7.1 used here, which has updated Ee ioni¬ 
sation temperatures). The fact that we find a peak at log T = 6.2 
is therefore not surprising. Our prominence is at a lower altitude, 
which means that there is more hot coronal emission along the 
line of sight than in the |Parenti & Vial| case. Also, the selection of 
EIS lines used means that our DEM is better constrained at this 
temperature than at lower temperatures, where there are limited 
lines available. The coolest line available for this DEM was O v 
at logT = 5.4. Below this temperature there are no constraints, 
and we see a sharp drop-off that therefore cannot be associated 
with any observational effects. 


We note here that in the range log T = 5.6 - 6.0 the slope of 
the tornado DEM (Eigure [T^ matches that of the 
DEM. 


Parenti & Vial 


7. Discussions and Conciusions 

This paper presents a number of plasma diagnostics using an 
EIS data set that captured a tornado above the solar limb. We 
find a similar line-of-sight velocity pattern to |Su et ar] ( |2Q14| ) at 
logT = 6.2. We show that this pattern in the line-of-sight ve¬ 
locities, suggestive of rotational motion, persists at temperatures 
down to logT = 6.0 with the Eeix 197.862 A line, and up to 
temperatures of logT = 6.3, with the Eexiii line at 202.044 A. 
Lines formed at temperatures lower than log T = 6.0 do not show 
any signs of a Doppler split about the tornado axis. It should be 
noted, however, that it may still be the case that this pattern ex¬ 
ists at lower plasma temperatures, but for reasons discussed in 
Section!^ we cannot see them with this data set. 
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Using the CHIANTI atomic database, we have been able to 
give an estimate of the electron density in the region of the tor¬ 
nado based on Fe xi and Fe xii lines, both of which are formed 
at around logT = 6.2. At this temperature we obtain an electron 
density of log = 8.5 when looking directly towards the centre 
of the tornado. In the surrounding plasma, however, we recover 
log He ~ 9. In other words, in all diagnostics available we see a 
dip in density when looking towards the tornado location, com¬ 
pared with the ambient corona. 

These results suggest that the hot moving plasma forms a 
sheath region of lower density around the cool core. On the edge 
of the tornado, we could be seeing denser, cooling threads. This 
picture seems consistent with the multi-threaded prominence 
model of |Luna et al.| ( |2012| . 

The non-thermal line widths inferred from our analysis may 
indicate broader line profiles for two Fe lines at the tornado lo¬ 
cation, though this result is not entirely certain. If this broad¬ 
ening is, however, localised to and caused by the tornado, then 
this could be explained by the presence of turbulence, or by the 
Kelvin-Helmholtz instability, or possibly other mechanisms. 

We provide the first prominence DEMs using FIS data. Pre¬ 
vious authors, using other instrume nts, found that prominence 
DEMs pea k at around logT = 6.0 ( |Wiik et al.|1993[ Parenti & 


Vial|2007) , whereas we find a peak at logT = 6.2. We attribute 


this difference to the selection of lines available from FIS, which 
are mostly formed at hot coronal temperatures, differences in 
the atomic data used, as well as differences in the relative alti¬ 
tudes of the prominences in question. We also find differences 
between the DEM measured in the tornado and that from the 
main prominence body. In the tornado there appears to be more 
contribution at lower temperatures (logT = 5.4) than we find in 
the prominence. We also find differences in the gradients of the 
tornado DEM, taken between logT = 5.65 and logT = 6.10, 
and the prominence body DEM, taken between log T = 5.88 and 
logT = 6.10. In the tornado we have a gradient of 2.6^^^, and in 
the prominence body it is 6.4^^-^. 

Following this work we aim to go on to investigate the mag¬ 
netic structure of these tornado features, as well as the plasma 
properties over a broader temperature range. Along with the 
present study, these should help build a more accurate picture 
of the nature of solar tornadoes and how they relate to the over¬ 
all prominence structure. 
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